Gamma-Rays from Dark Matter Mini-Spikes in M31 
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The existence of a population of wandering Intermediate Mass Black Holes (IMBHs) is a generic 
prediction of scenarios that seek to explain the formation of Supermassive Black Holes in terms of 
growth from massive seeds. The growth of IMBHs may lead to the formation of DM overdensities 
called "mini-spikes", recently proposed as ideal targets for indirect DM searches. Current ground- 
based gamma-ray experiments, however, cannot search for these objects due to their limited field of 
view, and it might be challenging to discriminate mini-spikes in the Milky Way from the many as- 
trophysical sources that GLAST is expected to observe. We show here that gamma-ray experiments 
can effectively search for IMBHs in the nearby Andromeda galaxy (also known as M31), where 
mini-spikes would appear as a distribution of point-sources, isotropically distributed in a ~ 3° circle 
around the galactic center. For a neutralino-like DM candidate with a mass m x = 150 GeV, up to 
20 sources would be detected with GLAST (at 5cr, in 2 months). With Air Cherenkov Telescopes 
such as MAGIC and VERITAS, up to 10 sources might be detected, provided that the mass of 
neutralino is in the TeV range or above. 
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I. INTRODUCTION 

The nature of Dark Matter (DM) is, more than 70 
years after its discovery, still an open problem. It is com- 
monly assumed that DM is made of Weakly Interact- 
ing Massive Particles (WIMPs), arisingin theories be- 
yond the Standard Model (see Refs. [H, for recent re- 
views), the most widely discussed DM candidates being 
the supersymmetric neutralino and the lightest Kaluza- 
Klein particle (LKP) in theories with Unified Extra- 
Dimensions [H, 0, Q. These particles will be actively 
searched for in upcoming high energy physics experi- 
ments such as the Large Hadron Collider (LHC, see e.g. 
Refs.[5|,la,LZ|,l8| for recent discussions in the context of DM 
searches), while hints on the nature of DM may already 
come from direct detection experiments aiming at detect- 
ing the nuclear recoils due to DM scattering off nuclei in 
large detectors (see e.g. @ for a recent update on the sta- 
tus of direct searches) . Alternatively, one could search for 
DM indirectly, i.e. through the detection of its annihila- 
tion products such as photons, neutrinos, positrons and 
antiprotons. The annihilation rate being proportional 
to the square of the DM density, ideal targets of indirect 
searches include all those regions where the DM density is 
strongly enhanced, due to gravitational clusterin g, a s in 
the case of the Galactic center jLu [lU [la Ua , Il4l. Ha. [1611 
and halo substructures [TJHlHSlHlillHlMtii, 
or because of energy losses capture in large celestial bod- 
ies, as in the case of the Sun and the Earth (see Ref. @ 
and references therein). 

Large DM overdensities can also form as a consequence 
of astrophysical pro cesses, such as the adiabatic growth 
of Supermassive |26l . l27l [28| or Intermediate Mass Black 
Holes [29|, [3(J. In fact, DM halos inevitably react to the 



'Electronic address: mfornasa@pd. infn.it, taoso@pd.infn.it, gianfranco. bcrtoneQpd. irifri.lt 



growth of black holes, leading, in the case of adiabatic 
growth, to the formation of large DM overdensities called 
spikes 26]. A DM cusp with a power-law density profile 
p ex r~ 7 , gets redistributed after the BH growth into a 
steeper profile p sp oc r _7sp , with j sp = (9 — 27)/ (4 — 7), 
within the radius of gravitational influence of the Black 
Hole (BH) (see below for further details). BHs can thus 
be thought as annihilation boosters, because the annihila- 
tion rate after their growth is boosted by several orders of 
magnitude, making these objects ideal tar gets for indirect 
DM searches. Even in absence of mergers [31(, and ignor- 
ing a possible off-center formation [32] , a spike around the 
Supermassive BH at the Galactic center would inevitably 
be destroyed by the combined effect of gravitational scat- 
tering off the observed stellar cusp at the GC, and DM 
annihilations [141 ] . The very same gravitational processes 
can still lead to the formation of moderate enhancements 
called crests (Collisionally REgenerated STructures), but 
these structures do not lead to significant enhancements 
of the annihilation signal [33| . Mini-spikes around In- 
termediate Mass Black Holes (IMBHs) are more promis- 
ing targets of indirect detection, since they would not be 
affected by these dynamical processes, and they should 
appear as bright point-like sources, which could be easily 
detected by large field of view gamma-ray experiments as 
GLAST [3J] and further studied with ground-based Air 
Cherenkov telescopes (ACTs) [H such as CANGAROO 
[H, HESS [U, MAGIC [13 and VERITAS 
Here, we further explore the mini-spikes scenario, and 
focus on the population of IMBHs in the Andromeda 
Galaxy (also known as M31), a spiral galaxy very sim- 
ilar to the Milky Way (MW), whose center is located 
784 kpc away from us. We compute gamma-ray fluxes 
from DM annihilations around IMBHs in M31, and show 
that the prospects for detection with GLAST are very 
promising: in an optimistic case (a neutralino with a 
mass m y = 150 GeV and annihilation cross section av — 
GLAST may detect up to 20 point- 
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like sources (at 5a and with a 2 months exposure), within 
3° from the center of Andromeda. The proposed observa- 
tional strategy appears particularly suited for ACTs like 
MAGIC and VERITAS, (M31 is in a region of the sky 
not accessible to HESS), since the main difficulty in the 
search for Galactic mini-spikes is that they cannot per- 
form deep full-sky searches, due to their limited field of 
view. In the case of mini-spikes in M31, ACTs can search 
for them by scanning a small region of w 3° around the 
center of M31, and an effective exposure of « 100 hours 
in this region would be sufficient to probe the proposed 
scenario, at least for DM mass in the TeV range. The 
next-generation Cherenkov Telescopes Array (CTA)[39[, 
is expected to significantly improve the sensitivity, in- 
crease the field of view and decrease the energy threshold 
with respect to existing ACTs, thus representing an ideal 
experiment for the proposed scenario. 

The paper is organized as follows: next section (Sec- 
tion [TlJ is devoted to describe the formation scenario of 
IMBHs. We then (Section [ff[J present the IMBHs cata- 
logue and how it is adapted to the Andromeda Galaxy. 
We compute the gamma-ray flux emitted by each point- 
like spike around an IMBH, considering a particular en- 
ergy annihilation spectrum for a DM particle. In section 
II V Al we estimate the prospects for detection for a generic 
ACT. Then in section HvTH we turn to GLAST. Finally 
our results are discussed in Section IVl 



II. INTERMEDIATE MASS BLACK HOLES 
A. IMBHs formation scenario 

IMBHs are compact objects with mass larger than w 
20M Q , the heaviest remnant of a stellar collapse and 
smaller than ss 10 6 M Q , the lower end of the mass range of 
SuperMassive Black Holes (SMBH) jtj. The theoretical 
and observational motivations for IMBHs were recently 
reviewed in Ref. (42J. For instance, Ultra-Luminous X- 
ray point sources (ULXs) could be interpreted as accret- 
ing IMBHs, since alternative explanations in terms of 
AGNs, neutron stars or SMBHs appear to be problem- 
atic or even ruled out [H, [43| . 

From a theoretical point of view, a population of mas- 
sive seed black holes could help to explain the origin 
of SMBHs. In fact, observations of quasars at redshift 
z w 6 in the Sloan Digital survey [H, [H, H(| suggest 
that SMBHs were already in place when the Universe 
was only ~ 1 Gyr old, a circumstance that can be un- 
derstood in terms of rapid growth starting from massive 
seeds (see e.g. Ref. [47]). 

In fact, a generic prediction of scenarios that seek to ex- 
plain the properties of the observed SuperMassive Black 
Hole population, is that a large number of "wandering" 
IMBHs should exist in DM halos 0, [H, Despite 
their theoretical interest, it is difficult to obtain con- 
clusive evidence for the existence of IMBHs. A viable 
detection strategy could be the search for gravitational 



waves produced in the mergers of the IMBH popula- 
tion [5l|, H2, [H, HH l55|. [561 . with space-based interfer- 
ometers such as LISA [57j |. 

In Ref. [3(j, two scenarios for IMBHs formation have 
been considered. The first posits IMBHs as remnants 
of the collapse of Population III stars. Zero-metallicity 
Population III stars are more massive than more recent 
metal-enriched stars and, if heavier than 260M Q , they 
would collapse directly into black holes (|42| and Refs. 
therein) . 

Here, we will focus only on the second scenario, based 
on Ref. [49j], where IMBHs form at high redshift from 
gas collapsing in mini-halos. If the latter are massive 
enough, proto-galactic disks form at the center of each 
halo. Gravitational instabilities introduce an effective 
viscosity that causes an inward mass and an outward 
angular momentum flow. The process goes on till it is 
interrupted by feedback from star formation (1-10 Myrs) 
that heats the disk. Then the so-formed object undergoes 
gravitational collapse into a black hole. A characteristic 
mass scale of 10 7 M Q is imprinted to the mini-halo by the 
requirements that it is heavy enough to form a gravita- 
tional unstable disc and that the black hole formation 
timescale is shorter than the typical major mergers one. 
The resulting black holes have a mass log-normally scat- 
tered, with a a, = 0.9, around the mean value of [43 ]: 




where n is that fraction of the baryonic mass which 
loses its angular moment that remains in the remnant 
black hole. / is the fraction of the total baryonic mass 
in the halo that has fallen into the disc, M V i r is the 
halo virial mass, z is the redshift of formation and t the 
timescale for the evolution of the first generation of stars. 

In Ref. [3(| Bertone et al. have studied the population 
of IMBHs in the MW, following the evolution of mini- 
halos hosting IMBHs at high redshift (populated with 
the prescriptions of Ref. [49||). down to redshift z = (see 
Ref. [3(| for further details). As a result, they obtained 
200 realizations of the IMBHs population in the Galaxy, 
that were used to produce 200 mock catalogs of DM mini- 
spikes, and to study the prospects for detection of these 
objects in the Galaxy. The average number of unmerged 
IMBHs was found to be N = 101 ± 22, and each of these 
objects is characterized by its mass, distance from the 
center of the galaxy, and surrounding DM distribution. 

B. DM distribution around IMBHs 

Following earlier work on the dynamics of stars and 
DM around compact objects (see Ref. [58[ and references 
therein) , Gondolo and Silk have shown that the adiabatic 
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growth of a massive black hole in the center of a dark halo 
modifies the distribution of the surrounding DM, induc- 
ing an enhancement of the density called "spike" [2(| . 
They focused their attention on the SMBH at the center 
of our Galaxy, but the same formalism can be applied 
also to IMBHs. The initial DM distribution in all mini- 
halos can be adequately parametrized with a Navarro, 
Frenk and White (NFW) profile 

«•■>-» (7) + < 2 > 

where r S) called the scale radius, sets the radius at which 
the profile slope changes. The new profile after the adi- 
abatic growth, will be (26l |: 

/ \ -7/3 

Ps P (r) = p(r sp ) ( — ) , (3) 

\ r sp / 

where p is the density function of the initial NFW profile. 
r sp gives the upper limit inside which Eq. [3] is considered 
valid and is related to the radius of gravitational influ- 
ence of the black hole r^: r sp « 0.2r/j [6(|, where is 
implicitly defined as: 

M(r<r h )= I p(r)r 2 dr = 2M. 
Jo 

with M, is the mass of the black hole. 

The spike profile diverges at low radii but annihila- 
tions set an upper limit to the physical density. Solving 
the evolution of the DM particles number density, one 
finds that the upper limit depends on the microphysical 
properties of the DM particles (mass and annihilation 
cross section) and on the evolution timescale of the black 
hole. We denote the distance where the density equals 
this upper limit rn m , and following Ref. [301 ] we define a 
cut-radius for our density profiles: 

r cut = Max[AR Schw ,r Hm ], (4) 

where Rschw = 2.95 km M/Mq is the BH Schwarzschild 
radius. The density between Rschw and r cu t is assumed 
to be constant to p sp {r cut ). 

III. GAMMA-RAYS FROM IMBHS IN M31 
A. IMBHs in M31 

Although similar, the Milky Way and Andromeda do 
not have exactly the same properties. The mock catalogs 
of IMBHs built for our Galaxy, thus have to be modified 
to account for the different average number and different 
spatial distribution in the host halo. A comparison be- 
tween the properties of Andromeda and of the Galaxy is 
shown in Table [I] 

We start from the mock catalogs obtained in Ref. [3(| 
and we rescale the total number of objects by the ratio 




Figure 1: Differential photon spectrum per annihilation. Dif- 
ferent parametrizations and annihilation channels are shown. 
Solid line (FPS) is an analytic fit relative to the bb channel, as 
obtained in Eq. [7] Dashed line (Eq. llO|l is relative to the same 
annihilation channel bb, but with a different parametrization 
of the FFs (see Eqs. H and [TO). Dotted line (BBEG) is rel- 
ative to B 1 annihilations and includes final state radiation 
from annihilation to charged leptons [651 ] (see text for more 
details) 



between the host halo masses, since the number of un- 
merged IMBHs scales linearly with the host halo mass, 
and the galactocentric distance by the ratio of virial radii. 
We obtain for M3I an average number of IMBHs per re- 
alization Nm3i = 65.2 ± 14.5. The mass spectrum re- 
mains unchanged, with an average mass around I0 5 Mq, 
while the average distance from the center of the galaxy 
is 32.31 kpc. We have verified that our rescaling proce- 
dure reproduces in a satisfactory way the properties of 
the IMBHs population in Andromeda, by comparing our 
results with a limited number of mock catalogs obtained 
as an exploratory study in Ref. 30] . 





Milky Way 


Andromeda 


Distance to the center [kpc] 
Virial Radius [kpc] 
Virial Mass [Mq] 
r a [kpc] 


8.5 

205 
1.0 x 10 12 

21.75 
5.376 x 10 6 


784.0 
180 
6.8 x 10 11 

8.18 
3.780 x 10 7 



Table I: Distance from the Sun (in kpc), virial radius (defined 
as the radius within which the density reaches 200 times the 
critical density, in kpc), virial mass (in solar masses) and the 
two NFW density profile parameters (in kpc and MQkpc -3 
resp ectively) , both for the MW and the Andromeda Galaxy 
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Figure 2: Left: Map of gamma-rays emission (above 100 GeV, in cm _2 s _1 ) by DM annihilation (m x = 1 TeV) around IMBHs 
in Andromeda, averaged over all realizations. Bins are 0.1° wide, to match the angular resolution of ACTs and GLAST. 
The circle shows for comparison the M31 scale radius r s . Right: Luminosity function of IMBHs (fluxes are in cm _2 s _1 ), for 
m x — 0.3, 0.5 and 1 TeV. The vertical line shows the contribution of the smooth component of the M31 halo, assuming a NFW 
profile and m x = 1 TeV. 



B. Gamma-rays flux from IMBHs in M31 

Once the mock catalogs of IMBHs in M31 have been 
obtained, it is possible to calculate the gamma-ray flux 
from each IMBH in every realization. The calculation 
follows Eq. 14 in Ref. 
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where $ = 2.7xl0 -14 cm^s" 1 , d is the IMBH distance 
to the observer, av is the DM annihilation cross section 
times relative velocity and m x is the DM particle mass 
(the letter x, usually adopted for neutralino, is used here 
to denote a generic WIMP candidate). r cut and r sp , rep- 
resent the inner and outer size of the spike, as discussed 
in the previous section. 

dNj (E) I dE is the differential photon yield per annihi- 
lation, that can be expressed as: 



dN 1 {E) 
dE 



a dN«{E) 



dE 



(6) 



where B a is the branching ratio BR(x% — > aa) and 



dNOj/dE the secondary photon spectrum relative to the 
annihilation channel ad. The latter term is thus a purely 
Standard Model calculation, while branching ratios have 
to be derived in the framework of new theories beyond 
the Standard Model, such as SUSY or UED. 

We review here different parametrizations of the pho- 
ton yield that have been recently proposed in literature. 
The first parametrization we focus on, has been obtained 
in Ref. [61( , and it is relative to annihilations into bb. The 
authors have parametrized the results obtained with the 
event generator PYTHIA [63[ as follows 



dN b Jx) 



dx 



^.a^b+cx+dx +ex L 



(J) 



where the parameters depend on the neutralino mass, 
and for the specific case m x — 1 TeV, (a, 6, c, d, e) = 
(-1.5,0.37,-16.05,18.01,-19.50). While for annihila- 
tion to rs 



dNZ(x) 

— f-- = x a {bx + cx 2 + dx 3 )e ex , 
dx 



(8) 



and for m x = 1 TeV, (a, b, c, d, e) = 
(-1.31, 6.94, -4.93, -0.51, -4.53). 

Alternatively, one can start from the most recent Frag- 
mentations Functions (FFs) (e.g. [H|]), describing the 
hadronization of partons into the particles of interest. 
The FF of b quarks hadronizing in neutral pions has been 
fitted with a simple analytic form that captures in a sat- 
isfactory way the behavior of the FF at large x finding 
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the following analytic fit 



/(*) 



7.53 



(9) 



Convolving the spectrum pions with their decay spec- 
trum into photons one finally obtains the differential pho- 
ton yield 



dN 7 (x) 
dx 



= I f(x')^dx'. 



(10) 



We have also considered an example inspired from the- 
ories with Unified Extra-Dimensions, where the role of 
DM is usually played by the first excitation of the hy- 
percharge gauge boson, and referred to as B^K Since 
the B^ annihilation into fermions does not suffer from 
chirality suppression, as in MSSM, we also include the 
contribution from annihilation to llj, as calculated in 
Ref. [65j, as well as the contribution from r fragmenta- 
tion, and usual from annihilations to bb, with the appro- 
priate branching ratio. The final state radiation arising 
from annihilation to charged leptons has a characteristic, 
very hard, spectral shape (65l . I6r3 | 



dN\{x 
dx 



l—e,fi 
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(1-x) 



(11) 



The three prescriptions for the annihilation spectrum 
are plotted in Fig. Q] (for m x — 1 TeV). As expected, all 
spectra are very similar up to x = E/m x ^0.1, but the 
spectrum relative to B^ annihilations is harder at large 
x, and exhibits a distinctive sharp cut-off at x = 1. 

To show the small effect that the adoption of differ- 
ent annihilation spectra has on the prospects for indirect 
detection, we have calculated the DM annihilation flux 
from the smooth component of the M31 halo, assuming 
a NFW profile with the parameters described in Tab. U 
above. The results are displayed in Table ITTI and as one 
can see, differences are within a factor of 2. In the re- 
main of this paper, we will thus work only with the first 
analytic fit, since the uncertainties associated with other 
astrophysical and particle physics parameters are signif- 
icantly larger. 

By calculating the gamma-ray flux in Eq.[5]for IMBHs 
in all realizations, we obtain a gamma-ray map, above 
Ethr = 100 GeV, from mini-spikes in Andromeda. Each 
M31 realization actually produces a different emission 





M31 flux [cm" 2 s" 1 ] 


FPS [61] 
Eq.[T0] 

BBEG [65] 


1.33 x KT 14 
9.79 x KT 15 
1.60 x KT 14 



Table II: Gamma-ray flux over 100 GeV from Andromeda 
(in cm~ 2 s _1 ) for a smooth NFW, and for the different 
parametrizations discussed in the text. Differences among 
the predicted fluxes are within a factor of 2. 



10-' 



io- ! 



*(<8) [cm" 2 s 1 ] 



icr" = 



10"" 



10" 1 



10" et 



10" 



signal 

Galactic/ExtraGal. back. 
Hadronic/Electronic back. 



en 



0.5 



1 1.5 



2.5 



Figure 3: Gamma-ray flux (in cm~ 2 s~ 1 ) from DM annihila- 
tion around IMBHs (solid thick line) , integrated over a cone of 
size 9 towards the center of M31, as a function of 9. We show 
for comparison the hadronic/electron background, assuming 
Eh — 10~ 2 (solid thinner line) and the diffuse extragalactic 
background (dashed line). 



map, so in Fig. [5] we show the average of all 200 maps, 
which clearly exhibits, as expected, a strong enhance- 
ment of the flux in the innermost regions of the Galaxy. 
The pixel size matches the angular resolution of ground- 
based telescopes such as VERITAS and MAGIC, and of 
GLAST. For the map, a DM particle mass of 1 TeV and 
an annihilation cross section av = 3 x 10 _26 cm 3 s^ 1 have 
been adopted. Note that this map alone does not provide 
any information on the detectability of the fluxes, which 
will be discussed in detail in the next 2 sections. Note 
also that the actual distribution of IMBHs, will provide 
(as we shall see later) a much more 'patchy' emission, far 
from the average smooth behaviour shown in Fig. [5] 

We also show in the right panel of Fig. [2] the lumi- 
nosity function of IMBHs (sum of all realizations), for 
different values of the DM particle mass. The distribu- 
tion is approximately gaussian, and the average flux of 
tMBHs is larger than emission due the smooth compo- 
nent. The dependence from the mass results in due to 



a balance between the 



-9/7 



dependence in Eq. [5j and 
the m x dependence of the upper limit in the integral of 
the energy spectrum. Having set in the figure an energy 
threshold E t hr = 100 GeV, the luminosity flux towards 
higher fluxes when the mass increase. We will come back 
later to this threshold effect, that leads to higher fluxes 
for higher masses when m x ~ E t h r despite the explicit 
m x dependence of the annihilation flux. Meanwhile 
we note that this effect disappears when m x 3> E t h r , as 
can be seen from Table IIII1 



6 




Figure 4: Left (right) panel shows a map of the gamma-ray flux in units of photons cm -2 s _1 , from DM annihilations around 
IMBHs in M31, relative to one random realization of IMBHs in M31. The size of the bins is 0.1° and the threshold for the left 
(right) panel is 100 GeV (4 GeV) as appropriate for ACTs (GLAST) . The circles highlight IMBHs within the reach of current 
ACTs for a 5a detection in 100 hours (within the reach of GLAST for a 5a detection in 2 months). 



IV. PROSPECTS FOR DETECTION 

As we shall see the prospects for detection de- 
pend on the expected or measured experimental per- 
formances, but also on the atmospheric and astrophys- 
ical backgrounds. We perform separate analysis for 
Air Cherenkov Telescopes and the upcoming gamma-ray 
satellite GLAST. 



A. Prospects for ACTs 

The calculations in this section are performed for 
a generic ACT, but they are particularly relevant for 
two specific experiments: MAGIC and VERITAS. As 
for HESS, being located in Namibia, it cannot detect 
gamma-rays from the direction of Andromeda. 

To determine the significance of the signal from an in- 
dividual mini-spike, as calculated in the previous section, 





Average flux [cm 2 s 1 ] 


m x = 50 GeV 
m x = 150 GeV 
m x = 300 GeV 
m x = 500 GeV 


5.26 x 10~ n 
7.65 x lO^ 11 
6.92 x 10~ n 
5.81 x 10 -11 



Table III: Average flux from IMBHs in all 200 realizations 



(in cm 

10- 26 CE 



), for different values of DM mass, av 
1 and E thr = 4 GeV. 



= 3 x 



we compare the number of signal photons, to the fluctu- 
ations of the background 



VT~- Aft 



jA eff (E,e)^dEd9 



d<S> bk 
dE 



(12) 



dEde 



where T is the exposure time, A e ff the effective area, 
AJ7 the solid angle, d^bk/dE is the total background 
differential flux. 

For Air Cherenkov Telescopes, the main background 
is due to hadrons interacting with the atmosphere and 
producing electromagnetic showers. Following Ref. [67j 
[2J|, we consider 



d* h f E \ 

dQ.dE ' X VGeV/ 



-2.74 



P 



cm 2 s GeV sr 



(13) 



The ratio of the number of hadrons misinterpreted 
as gamma-rays, over the total number of cosmic ray 
hadrons, e^, provides an estimate of the telescope po- 
tential to discriminate the gamma-ray signal from the 
hadronic background. We adopt a typical value e/, = 
10 -2 , following Refs. [24| EjM]. The electronic contribu- 
tion to the background is [24j |: 



d<P e , / E \ 



-3.3 



cm 2 s GeV sr 



(14) 



and it is typically sub-dominant at the energies of interest. 

In Figure [3] we compare the DM annihilation signal 
with the different sources of background, as a function of 
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Figure 5: Number of detectable mini-spikes in M31 with GLAST (2 months) and ACTs (100 hours) as a function of the DM 
particle mass (left) and as a function of the angular distance from the center of M31 (right). In the left panel, error bars denote 
the 1 — a scatter among different realizations. In the right panel, the total number of objects is shown as an empty histogram, 
while the vertical lines denote the size of the region that contains 90% of the detectable IMBHs. 



the field of view. The minimum flux for a 5a detection 
with an effective area of A e f / = 3 x 10 4 m 2 [33] and an ex- 
posure time of 100 hours, is <f) m in = 1-6 x 10 _12 cm _2 s _1 . 

To produce this estimate we have considered values of 
effective area and angular resolution similar to MAGIC 
and the result is consistent with earlier estimates of 
the MAGIC sensitivity 0. An actual estimate of 
the instrument performance suggests that the mini- 
mum flux can be up to an order of magnitude higher [6!| . 

The fluxes in Fig.[^]are thus found to lie below the min- 
imum detectable flux so determined, so one may naively 
conclude that there is no hope to detect them with this 
experimental setting. However, for a more careful assess- 
ment of the prospects for detection, we need to estimate 
the detectability in each realization and then average over 
all the realization, and not viceversa. In fact, as demon- 
strated in the left panel of Fig. [H in each random real- 
ization, the emission is much more patchy, with a large 
number of high emission peaks, corresponding to indi- 
vidual mini-spikes, that can thus be resolved with the 
adopted angular resolution. Black circles highlight the 
position of objects brighter than the experimental sensi- 
tivity (indicated in the color scale by the black line). In 
total, for m x = 1 TeV, the number of detectable IMBHs 
is N^ a = 5.2 ±3.1, where the error is relative to the l-u 
scatter among different realizations. 

We note that current simulations indicate that the 
next-generation Cherenkov Telescopes Array (CTA)(39|. 
may significantly improve the sensitivity, down to <j) m in ~ 
10^ 13 cm _2 s~ 1 , thus leading to a substantial improve- 



ment in the prospects for detection. 



B. Prospects for GLAST 

The space satellite GLAST is expected to play a crucial 
role in indirect DM searches, thanks both to its ability 
to perform observations at energy scales comparable to 
the mass of common DM candidates and to its poten- 
tial of making deep full-sky maps in gamma-rays, thanks 
to its large (~ 2.4 sr) field of view (34[. Despite the 
smaller effective area, it is not affected, being a satellite, 
by the atmospheric hadronic and electron background. 
Furthermore, its lower energy threshold (30 MeV) allows 
to probe lighter DM particles, typically leading to higher 
fluxes. The angular resolution of GLAST is w 3° in the 
energy range 30 MeV-500 MeV, becomes 0.5° from 500 
MeV to 4 GeV, and reaches 0.15° above 4 GeV (H- 

As in the case of ACTs, we compare the expected fluxes 
with the photon background, which in this case, since 
GLAST will perform observations above the atmosphere, 
is mainly due to diffuse gamma-ray emission. The galac- 
tic and extragalactic background has been measured in 
0,[Dl by EGRET in the energy range between 30 MeV 
and 10 GeV and we extrapolate it to higher energies by 
fitting with a power-law with spectral index of -2.1. The 
resulting formula is 

d^extra/gal = 6 / ^J]'^ - 

dfldE ' \GeV J cm 2 s GeVsr' 

(15) 
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We note here that a large fraction of the observed 
gamma-ray background might be actually due to DM 
annihilations [72, [zl, [3 [ill, in particular if as- 
trophysical processes can boost the annihilation sig- 
nal [Zg, [ZZl- ln this case, the smoking-gun for this 
scenario would be the distinctive shape of the angular 
power-spectrum of the background, that may allow, al- 
ready with GLAST, the discrimination against ordinary 
astrophysical sources [zl| . 

The sensitivity above 30 MeV, i.e. the minimum de- 
tectable flux for a 5er detection with an exposure of 2 
months, is found to be <f> m in = 3.2 x 10 _8 cm _2 s _1 . This 
value is derived from Eq. [T2l adopting values of the en- 
ergy dependent effective area provided by and is 
consistent with GLAST sensitivity maps obtained in Ref. 
[80| . The integral flux above threshold from IMBHs, av- 
eraged among realizations and integrated in a 3° cone 
towards M31, is 3O = 1.3 x lO^cm^s" 1 . 

In the right panel of Fig. HJ we show the results of our 
analysis relative to a random realization, and adopting 
m x = 150 GeV, and av = 3x 10~ 26 cm 3 s _1 . Mini-spikes 
appear as high emission peaks, and can be easily resolved 
by selecting photons above 4 GeV, so that the angular 
resolution of GLAST approaches 0.1 degrees. Black cir- 
cles highlight those objects that produce a flux detectable 
at 5cr with GLAST, with a 2 months exposure. 

V. DISCUSSION AND CONCLUSIONS 

Although we have performed the analysis of the 
prospects for detection with GLAST and ACTs for 2 dif- 
ferent benchmark scenarios (essentially high DM particle 
mass for ACTs, low m x for GLAST), the analysis can be 
easily extended to any value of the particle physics pa- 
rameters of the annihilating DM particle. To explore the 
dependence on m x , we show in the left panel of Fig. [3] the 
number of objects that can be detected with the afore- 
mentioned experiments, as a function of the DM parti- 
cle mass. Near the experiment threshold, fluxes increase 
with mass. When m x ^> E t hr this threshold effect dis- 
appears and one recovers the expected behavior (smaller 
fluxes for higher masses). 

Similarly, one can plot the number of detectable ob- 
jects as a function of the angular distance from the cen- 



ter of M31, to estimate the region where most mini-spikes 
can be found. This is shown in the right panel of Fig. O 
where the total number of objects is also shown for com- 
parison. Vertical lines denote the angular size of the re- 
gion that contains 90% of the detectable IMBHs for the 
various experiments, which has a characteristic size of 
9 = 3.3°. 

We stress that while in the case of Galactic IMBHs 
the identification of mini-spikes will require a case-by- 
case analysis of the spectral properties of unidentified 
gamma-ray sources, the detection of a cluster of sources 
around the center of Andromeda would per se provide 
unmistakable evidence for the proposesd scenario. 

In conclusion, we have computed gamma-ray fluxes 
from DM annihilations in mini-spikes around IMBHs in 
the Andromeda Galaxy. We have studied the prospects 
for detection with Air Cherenkov telescopes like MAGIC 
and VERITAS and with the GLAST satellite, and found 
that a handful of sources might be within the reach of 
current ACTs, while the prospects for the planned CTA 
are more encouraging. The obvious advantage of the pro- 
posed scenario with respect to mini-spikes in the MW, is 
that they are not randomly distributed over the sky, but 
they are contained, at 90%, within 3 degrees from the 
center of Andromeda, and can thus be searched for with 
ACTs by performing a deep scan of this small region. 

The prospects for GLAST appear more promising, 
since an exposure time of 2 months allows the detection 
of up to of w 20 mini-spikes, that would be resolved as 
a cluster of point-sources with identical spectra, within 
a ~ 3° region around the center of Andromeda. Such 
a distinctive prediction cannot be mimicked by ordinary 
astrophysical sources. As in the case of IMBHs in the 
MW, null searches would place very strong constraints 
on the proposed scenario in a wide portion of the DM 
parameter space. 
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